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Abstract. The anomaly-induced effective action is a useful tool for deriving the 
contributions coming from quantum effects of massless conformal fields. It is well- 
known that such corrections in the higher derivative vacuum sector of the gravitational 
action provide the same exponential inflation (Starobinsky model) as the cosmological 
constant term. At the same time, the presence of a classical electromagnetic field 
breaks down the exponential solution. In this paper we explore the role of the anomaly- 
induced term in the radiation sector and, furthermore, derive the "equation of state" 
and the scaling laws for all terms in the Einstein equations. As one could expect, the 
^- ■ scaling law for the vacuum anomaly-induced effective action is the same as for the 

cosmological constant. 

m ' 
in 

O ' 1 Introduction 

O' 

It is well-known that the conformal anomaly is useful for various applications Quantum Field 
Theory. In particular, the anomaly-induced effective action has been explored in the cosmological 
setting about three decades ago pQ. Soon it was discovered that, in the absence of matter fields or 
radiation, the quantum anomaly-induced contributions lead to the Starobinsky model of inflation 
[2] (see also HJ [5J EJ El [8] for an alternative work and further developments). The traditional 
version of this inflationary model is based on the unstable exponential solution [2], that implies 
some special choice of the number of the quantum fields with different spins [7j. For instance, 
the present-day Universe with (presumably) only photon being an active quantum field, or the 
early universe where the active quantum content is described by the Minimal Standard Model of 
particle physics, satisfy the condition of unstable inflation. 

An alternative possibility is to consider the supersymmetric matter content of active fields, 
that leads to the stable inflationary solution at the initial stage of inflation. The transition 
from stable to unstable inflationary regimes can be associated to the decoupling of the massive 
s-particles [91 [10] in the universe where the inflation is slowing down because of the quantum 
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effects of massive fields [TO]- The same effect holds in the presence of the cosmological constant, 
which actually plays only a small role in this story [11| [T2T I13j . In both cases of stable and 
unstable inflationary solutions one usually assume that the universe is empty, that means there 
are no matter fields and/or radiation. 

All the mentioned massive or massless fields of different spins are virtual ones, they manifest 
themselves only through their contributions to the vacuum action. If the real radiation is present, 
there is no exponential solution for the conformal factor a{t) and the last tends to the correspond- 
ing FRW solution if the particle content corresponds to the unstable case and if the initial data 
are chosen in an appropriate way [I] . After a while, the dynamical system describing the universe 
enters the regime where the effect of higher derivative terms becomes negligible [6j [11], and 
the behaviour of the conformal factor is essentially the same as in the purely classical universe 
dominated by classical radiation content. 

The above conclusion is based on the analysis of the theory with the action which includes 
Einstein-Hilbert term, classical radiation and higher-derivative anomaly-induced gravitational 
contributions. However, there may be one missed component in this consideration. In the case 
when the background radiation is present, one may need to take into account also the anomaly- 
induced contribution to the electromagnetic part. Indeed, the classical radiation does decouple 
from the equation for the conformal factor, and the radiation density manifests itself only via 
the first integral of this equation, that is the first Friedmann equation. So, it looks interesting to 
check what is the effect of the anomaly-induced term, e.g., whether it is capable to produce some 
significant change in the acceleration of the universe. The same problem can be explored, also, 
for the radiation-dominated epoch afterthe inflation ends. 

One can consider a bit more general formulation of the problem. We can derive the "equation 
of state" for all components of the gravitational action, namely the Einstein-Hilbert term, cos- 
mological constant term, free radiation, quantum anomaly-induced contribution to the radiation 
part and quantum anomaly-induced contribution to the vacuum part. Then we can check how 
the corresponding "energy densities" depend on the scale factor. For instance, the comparison 
of these dependencies for the anomaly-induced vacuum terms and the cosmological constant can 
better explain why the two kind of vacuum actions produce similar exponential behaviour. This 
issue may be also interesting in view of the recent attempts to deal with the cosmological constant 
problem by taking the anomaly-induced contributions into account |14t [T5] . 

The paper is organized as follows. In Sect. 2 we write down the anomaly- induced terms in 
both gravitational and electromagnetic sectors and consider the relations between timelike and 
spacelike components of the diagonalized equations for the metric (which are Energy-Momentum 
Tensors in the electromagnetic field case) . These relations can be seen as equation of state for all 
the terms in the modified Einstein equations in the cosmological setting. In Sect. 3 we explore 
what is the effect of the anomaly-induced electromagnetic term for the rate of expansion of the 
universe in the two different situations, namely when the higher derivative metric dependent 
terms are present or not. Finally, in Sect. 4 we draw our conclusions and discuss some possible 
applications of the results. 



2 Classical and anomalous terms in the effective action 

The conformal anomaly is the typical theoretical phenomenon for massless conformal invariant 
quantum fields on some nontrivial external background. In case of massless conformal fields 
the action of vacuum (gravitational one) has to include, at least, the conformal invariant higher 
derivative part (see, e.g., [TB] for the introduction and [17J for a recent review of Quantum Field 
Theory in curved space) 

S HD = I d 4 xV=g(aiC 2 + a 2 E + a 3 V 2 R). (1) 

Here C 2 = & p-2Ei*0 + (1/3) R 2 is the square of the Weyl tensor and E = R* ua p - 4R 2 a/3 + R 2 
is the integrand of the Gauss-Bonnet topological term. The terms in the Lagrangian of JT]) satisfy 
the conformal Noether identity and, furthermore, do not affect the dynamical equation for the 
conformal factor of the metric. At the quantum level, however, the conformal symmetry is violated 
and this also affects the cosmological solution. 

In the cosmological setting, massless conformal invariant quantum fields corresponds to the 
early epoch when the energy of the photons is much greater than masses of at least some of 
the charged spinor fields. This condition can be easily satisfied in the inflationary period, es- 
pecially in the framework of the Starobinsky model, which has, usually, very high values of the 
typical energies at the end of inflationary period. Furthermore, this condition can be fulfilled in 
the radiation-dominated period after inflation, where many massive fields approximately can be 
approximately treated as massless. 

2.1 Anomaly-induced terms 

Consider the approximation of massless fields. In case of both gravitational and electromagnetic 
background fields, the conformal anomaly has the form 

< T£ > = - (wC 2 + bE + cV 2 R + /3F 2 ) , (2) 

where F 2 = F 2 is square of the strength tensor of the electromagnetic fields, w, b, c are the 
/3-functions for the parameters of the vacuum action and j3 is proportional to the electromagnetic 
charge /3-function. At one loop order, using the Minimal Subtraction scheme of renormalization, 
we get 

e - Sr^ N > " «=??"• (3) 



I 



as a sum over charged fermions and scalars with the multiplicities Nf and N s correspondingly. 
The one- loop values of w,b and c can be found, e.g. in [T8 j [TB I [TT] . 

It is well known that taking into account the conformal anomaly in the cosmological case leads 
to the Starobinsky exponential solution [2] for the conformal factor, if there are no matter fields. 



The detailed discussion of the anomaly- induced action of electromagnetic field and its relation to the more 
general result coming from the physical renormalization scheme can be found in [19] , see also [201 121| . 



At the same time, if the radiation is present, there is no such solution. One can naturally ask 
whether the anomalous electromagnetic term in ([2]) can change this situation. And more general, 
whether this term can affect the expansion of the universe at the early stage of its history. 

In order to use field quantities in the cosmological setting, one has to perform some space 
averaging. Obviously, {F 2 ) ~ {E 2 ) — {H 2 ) equals zero for a free radiation. But this does 
not apply, e.g., to the radiation-dominated early universe, because in this case there is also a hot 
plasma of other particles and the content of the universe does not reduce to a free electromagnetic 
radiation. As a qualitative simplest estimate we shall suppose that F 2 ^ and set its scale-factor 
dependence in accordance to its conformal property. One can assume, for instance, that at some 
fixed scale the magnitude of this term is proportional to the p®, that is the classical radiation 
energy density. This radiation density is supposed to describe not only electromagnetic fields, 
but also a hot plasma which fills the Universe. It is important to note that such non-trivial 
material content of the universe is indeed possible at the last stage of the stable inflation, where 
we observe oscillations of the conformal factor [10] and, consequently, production of photons and 
charged particles. 

The anomaly-induced effective action can be easily derived as a functional of the new variables 
g^ v and a, where g^ v = g^ u ■ e 2a and the metric g^ has fixed determinant. Disregarding the 
conformal invariant term in the effective action we arrive at the following expression |22| : 



f = f d 4 x^{waC 2 + ba{E--V 2 R) + 2baAa+ f3oF 2 } - 3c + 2b [d 4 x^/=g~R 2 , (4) 
J 3 36 J 

where F 2 = g^ a g u ^ F a pF a p = e~ 4cr F 2 . The expression (jl|) is the quantum correction to the 
classical action of vacuum. Let us note that the covariant (nonlocal and local) forms of the 
anomaly- induced action are well-known [221 [19] [20] , but the eq. (j4]) is sufficient for our present 
purposes. The total action has the form 

St = Seh + T + conf . invariant terms , (5) 

where the conformal invariant terms include the classical actions of radiation and of hot charged 
particles and Seh is the Einstein-Hilbert term 

Seh = - j^ j d 4 x^g~ (R + 2A) . (6) 

In the case of the cosmological, FRW metric, the classical massless fields decouple from gravity. 
Taking into account the conformal properties, it proves useful to rewrite the expression §5§ in a 
more detailed form 

St = Seh + Shd + S r + Thd + Tp . (7) 

Here Shd and S® are classical higher derivative metric and radiation (including massless charged 
fields) conformal invariant actions. Thd and Tn are parts of anomalous action (jl|). In what 
follows we will use the same indications for all quantities, including the trace of the stress-tensor, 
Ti = (T E H,T HD ,f HD ,T^,Ti3), energy density p { and pressure p i: with T i = p i - < i-p i in the 
corresponding reference frame. On the top of that, we will sometimes use notation for the total 
expression in the radiation sector, like p r = p° + pp. 



2.2 Energy density and pressure 

Consider the stress-energy tensor, whose components are given by the variational derivative of 
the total effective action (J7J), 

T# = - i_ * St . (8) 

V~9 d 9a(3 

In order to calculate separately the contributions of the Einstein-Hilbert term from the terms 
with high derivatives (HD) which come from the quantum contributions and the electromagnetic 
one, we present the trace of the stress-energy tensor in the form 

1 A Q 

T = -* -p- = T E h + if + T HD + f HD + Tp = p t - 3p t , (9) 

or da 

where 

Pt = PEH + Pr + PHD + PHD + P/3 = ~ 1q (10) 

and (in Cartesian coordinates) 

Pt = PEH +P°r+ PHD + PHD +Pp= T? = T| = T* (11) 

are density-like and pressure-like components, correspondingly. The indices and bars of all quan- 
tities are in accordance with the ones of the actions in the r.h.s. of ([7]). We find useful to introduce 
such notations even for the Einstein-Hilbert term, despite the physical sense of the quantities is 
different in this case (as it is, of course, for the higher derivative vacuum terms, too). 

Our purpose is to see how all p's and p's depend on the scale factor a(t) and also how they 
behave during and after the inflationary period. One can find the densities for the components 
by assuming that the conservation law is satisfied separately for each of the stress-energy tensors 
Ti = (Teh,T®,Thd,Thd,T/3) in ([9]). In terms of the cosmic scale factor the conservation law 
can be expressed as 

rri 

d (pi a 3 ) = -pi d(a 3 ) , where pi = -^— ^ • ( 12 ) 

Following this standard procedure, we can immediately see that since the trace of the classical 
radiation stress tensor is zero, the equation of state is as it is supposed to be, 

T r ° = 3p° r - p° r = , hence p° r = ^ . (13) 

Correspondingly, this term does not contribute to the equation of motion for a(t). 

The first observation is that, since the trace is zero for the Shd term, the corresponding 
equation of state is exactly the same as for the free radiation, phd = Phd/3- F° r other three 
terms we obtain, in terms of conformal time r] (as usual, dt = a(r/)drj), 

m 3 \a" 2A1 , , 

Te « = ^g y-T ■ (14) 
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J3F 2 



(15) 



(16) 



Using Eqs. (|12|) . the solution for all the densities pi come from the differential equations of 
the form 



dpi _ ^fH_ _ Tj_ 
da 3 3 a 3 3a 3 

A general solution for this non-homogeneous equation (|17p is 

Pi (a) = C(a)a~ A , 

where the coefficient C{a) is obtained by the integration of 

dC T 3 / 
— = Tj a a . 
dry 



(17) 



(18) 



(19) 



Integrating (|19j) for each of remaining stress-energy tensor components above and substituting 
them into Eq. (|18p . we arrive at the following results: 
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Indeed, the equations (|2ip derived here are well-known, they are exactly the same as the ones 
obtained in pQ , and also recalculated in [8] . It is easy to see that these formulas are quite different 
from the ones for the cosmological constant in (|20p . One can expect that for the general form of 
0(77) these two different equations of state will definitely produce different contributions. However, 
in short we will see the effect of the two terms is equal for the exponential inflation case. 



3 Cosmological solutions with anomalous terms 

Here we consider the effect of radiation anomalous term on the behaviour of the conformal 
factor of the metric, in the framework of a FLRW cosmology, with and without higher derivative 
anomalous terms. 
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3.1 Stable anomaly- induced inflation with radiation term 

As a first step, consider the equation of motion including the higher derivative terms. The most 
useful choice of variable is the conformal factor as a function of cosmic time, cr(t). The last 
quantity is defined as a = In a. The equation of motion can be obtained from the 00-component 
[H [2] or directly from the trace T = p t — 3pt = [7]. In terms of r = t/tpi = Mpit, where tpi is 
the Planck time unit tpi ~ 5.3 x 10 -44 s, the equation has the form 

/ b 
a + 7 aa +4 3 

V c 

~ ^(-+2- 2 -^)--(^4-)e— = 0. (23) 

In this equation the contribution of the cosmological constant term is written in terms of vacuum 
energy density p A = A/(8vrG) = AM pl , where M Pi = l/y/8irG = 2.44 x 10 18 GeV is the reduced 
Planck mass. 

The direct inspection shows that, in the presence of the /3F 2 -term, there is no exponential 
solution. This demonstrates that the importance of vacuum for such solution hold also when we 
take the anomaly in the radiation sector into account. 

The next question is what is the role of the /3F 2 -term for the case of a stable inflation. 
In principle, one can expect two different situations: (i) The anomalous term slows down the 
exponential inflation, as it happens with the terms generated by the quantum effects of massive 
light fields |10t lll|; (ii) The anomalous /3-F 2 -term decreases very fast and soon becomes negligible. 
The numerical analysis show that this last behaviour actually takes place. For the illustration we 
present the corresponding plots for the case of Minimal Supersymmetric Standard Model (MSSM) 
in Fig. 1. As we have already mentioned in the Introduction, the supersymmetric particle content 
is the most interesting here, because it provides stable inflation, making the possible effect of the 
radiation term (or the absence of such effect) the most explicit and clear. 

In order to find the cosmological evolution of all the densities, we solve eq. (I23p numerically 
and then replace the solution into the expressions which directly follow, in particular, from (|20p . 
(pZTj) and (f22j) . We present these results, in Fig. 1, as functions of r = t/tpi, where 

-^r =3a ~wr (24) 

P %r)=e-^p r (r = 0), (25) 

PHd(t) ( .... .2- 1 -2 b .A 

—^- = -6c(aa+3oa--o -- a j , (26) 



Mi, M\ 



ae" 4 * 7 , (27) 
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for the MSSM particle content with N 1A / 2fi = (12,48, 104). 
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Figure 1: We have assumed here the MSSM particle content with -/V^^o = (12,48,104) and 
took the numerical value f3F 2 = —0.1. The plots for pud and peh rapidly tend to constants, 
because the inflation is stable and in the exponential regime the curvature components behave like 
constants. In this regime, it shows exactly the same scaling law as the density of the cosmological 
constant. 



3.2 Radiation-dominated evolution after inflation 

In this section we shall investigate the effect of quantum corrections at the period after the 
anomaly- induced inflation ends and the higher derivative terms in Q become negligible. Then 
the relevant part of the total action has the form 

/ **^ R + S — ~ + / *V=I/^ , (28) 



16ttG 

where S c i ass . ma tter is the classical action of the matter fields. We are interested in the period when 
matter and radiation are very hot and can be treated as conformal. Then, the classical massless 
fields in S c i aS s.matter decouple from the conformal factor of the metric and the effective equation 
of motion, in terms of conformal time, has the form 

A a " + ^ = 0. (29) 



4ttG 



Let us, as before, denote the derivative with respect to the cosmic time t by a point. Then 
Eq. (f29j) becomes 

• 2 9 

„ + --£,=<>. (30) 



In the last equation we have introduced a useful notatiorQ 

2 8ttG(3F 2 
V> = 3 • ( 31 ) 

As a first step in solving Eq. ()30p we find the relation between a(t) and H = a/a 



H2 = C ±J Mna^ (32) 



where C is an integration constant. In order to clarify the sense of this constant, let us consider the 
standard classical model with /i = 0. In this case, from Eq. (|32p follows a(t) = C(t — to) 1 , where 
to is some fixed instant of time. Different choices of to can be compensated by the renormalization 
of C, so we set to = 0. On the other hand, by solving the Friedmann equation we obtain 

r327rGp°(t = 0)iV4 r 

a{t)=[ ^y -J -a -Vi, (33) 

where p®(t = 0) and ao are the energy density of the electromagnetic field and the scale factor 
of the metric at the instant t = 0. The comparison of the two expressions for a{t) lead to the 
relation C = Pr } — ' . It is natural to fix that the ao corresponds to the fiducial metric g^ u . In 
what follows we put ao = 1. Then the elements of the solution (f32l) satisfy the relation fi 2 <^ C 
if the anomalous contribution and energy density satisfy the relation \(3F 2 \ <C 4p}. . As we have 
indicated above, this relation is quite natural, because for the free radiation F 2 = and the 
presence of the anomalous term is due to the interaction with other fields which have energy 
density much smaller than the one of radiation. 



'Let us note that F 2 = 2({H 2 ) — {E 2 }) is typically positive and that, according to our definition @, /3 < 0. 



Finally, the general analytic solution of (|30p can be presented in the form 

t = — / e 2z dz = J- Erfi(V2a + 2C/^ 2 ) - Erfi(V2C/ju) . (34) 

fj, J p y z L 

The disadvantage of this formula is that it becomes singular in the classical limit p, — >■ 0. In order 
to solve this difficulty, one can derive an aproximate solution by treating the term with p as a 
small perturbation, 

VC^f(l-^-lno). (35) 

As one can see from the last relation, the expansion of the universe performs slightly faster as a 
result of the quantum effects related to the electromagnetic anomalous term f3F 2 . 

When the universe expands, the radiation temperature is decreasing. It is instructive to 
find the temperature relations when the quantum term in the solution (|35[) is relevant. The 
lower bound for the relevant temperature is defined by the energy corresponding to the moment 
when the lightest charged fermion decouples and the upper bound is the energy scale when the 
inflation ends and the higher derivative terms in Eq. ([4]) become negligible. In the framework of 
the modified Starobinsky model [91 \TU[ [TT] , the scale of the graceful exit from the anomaly- induced 
inflation depends on the scale of the supersymmetry breaking, and may vary from H = 300 GeV 
to H = 10 14 GeV for different gauge theories. Let us notice that, contrary to the case of a purely 
gravitational background, the mixed electromagnetic-gravitational background corresponds to the 
loop diagrams with external lines of both electromagnetic potential and metric perturbations. 
Indeed, the decoupling scale is defined by the energy of the photons which is much greater than 
the energy of the gravitons. 

In order to find the temperature of the radiation corresponding to f)35f) . we can use the 
Friedmann equations 

a 2 , 8vrG 2a a 2 

- = H 2 = —p r , — + - = -87rG Pr , 36 

a z d a a z 

where p r is the radiation pressure and the thermodynamical relation is close to p r ~ p® = jg T 4 . 
Then we arrive at the formula 

T 4 = ™ H 2 ■ M 2 Pl . (37) 

The next requires the equation of state for the anomalous term, which was obtained in the 
previous section. Using this result directly leads us to 



C + p 2 lna 

8ttG o 1 



Pr : = ^^7 4 38 



and 



_ C + p 2 (lna-l) 1 1\P\F 2 

Pr ~ SvrGa 4 " 3 Pr + 3 a 4 " ^ J 
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As one can see here, the quantum effect on the background electromagnetic fields decreases the 
radiation pressure. The dependence between the temperature and the scale factor is given by 



a 



45(C + /i 2 lna 
8^ 3 G 



1/4 1 r 15 



a 



1/4 



- (p° r - \/3\F 2 lna) . (40) 



Naturally, the quantum effects produce some deviation from the usual classical formulas, namely, 
the p r is a bit larger than p® for a given temperature. 



4 Conclusions and discussions 

We have considered a cosmological applications of the vacuum quantum effects in the radiation- 
dominated universe and found that the J3F 2 term in the conformal anomaly leads to a slight 
modification of the evolution law and the thermal history of the universe. In the transitional 
period between inflation and radiation dominated universe the f3F 2 gives a non-zero contribution 
to the acceleration of the universe, that is different from the classical radiation. It would be 
interesting to explore further physical consequences of this effect. 

The relation for the anomaly-induced effective action in the radiation sector can be useful 
for investigating the general features of the gravity with anomaly- induced quantum corrections. 
In particular, it would be very interesting to explore the stability of the corresponding semi- 
classical solution, for a realistic particle content, that means non-stable Starobinsky inflation. 
This problem is well-known as a problem of stability of Minkowski and de Sitter spaces (see, e.g., 
|23[ [24~1 [25] and further reference therein) . Our previous analysis also shows that the stability 
conditions for the conformal factor fo the metric may be different for the flat space from one side 
and for the dS space from another one [11 j . It would be very interesting to check out what are 
the conditions of stability of the classical solution in a general case, for different stages of the 
universe expansion. The effect of radiation in the anomaly-induced action (refgeneral solution) 
is potentially relevant on this respect. We hope to report on this issue in a clase future. 
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